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Abstract —In this paper, we consider a mixed RF and hybrid 
RF/FSO system where several mobile users transmit their data 
over an RF link to a relay node (e.g. a small cell base station) and 
the relay forwards the information to a destination (e.g. a macro 
cell base station) over a hybrid RF/FSO backhaul link. The relay 
and the destination employ multiple antennas for transmission 
and reception over the RF links while each mobile user has a 
single antenna. The RF links are full-duplex with respect to the 
FSO link and half-duplex with respect to each other, l.e., either 
the user-relay RF link or the relay-destination RF link is active. 
For this communication setup, we derive the optimal resource 
allocation policy for sharing the RF bandwidth resource between 
the RF links. Our numerical results show the effectiveness of 
the proposed communication architecture and resource allocation 
policy, and their superiority compared to existing schemes which 
employ only one type of backhaul link. 

I. Introduction 

Free space optical (FSO) has recently re-emerged as an 
attractive option for wireless backhauling due to its large 
usable bandwidth compared to traditional radio frequency 
(RF) backhauling m. Moreover, FSO transceivers employ 
very narrow beams which leads to several interesting features 
such as secure and interference-free communication. However, 
these beneficial properties of FSO systems come at the expense 
of some drawbacks and challenges. One of the main challenges 
is the unpredictable connectivity of the FSO link due to 
atmospheric impairments such as atmospheric turbulence and 
visibility limiting conditions including snow, fog, and dust ID. 

A possible strategy to mitigate the unpredictable connectiv¬ 
ity of the FSO link are so-called hybrid RF/FSO links where 
an RF link is employed as a back-up for the FSO link ||2l, 
0. This option is attractive particularly due to the fact that 
distortions in the RF and FSO links are caused by different 
phenomena. Hence, a hybrid RF/FSO system is much more 
likely to maintain connectivity than a pure FSO system. FSO 
transceivers are mainly attractive for nodes with fixed locations 
where a line-of-sight link can be established. Motivated by 
this limitation of FSO for networks with mobile nodes, mixed 
RF/FSO systems were proposed in the literature. Here, the RF 
and FSO links are cascaded, i.e., the mobile nodes employ 
RF links to send their information to an intermediate fixed 
node, i.e., a relay node, and the intermediate node forwards the 
information to the final destination via an FSO backhaul link 
a-ia. This communication setup can model several practical 
applications including: i) Cellular communication where the 
mobile nodes send their data to a relay station and the relay 
station forwards the data to the base station; ii) Femto cell 
networks where the mobile nodes in a building floor send 
their data to an access point and the access point forwards 
the information to the macro base station. The performance 
of mixed RF/FSO systems was investigated in a, a for a 
single-user amplify-and-forward (AF) relay network, and in 
0 for a multi-user decode-and-forward (DF) relay network. 


In general, the end-to-end performance of dual-hop com¬ 
munication is limited by the weakest link. Hence, in mixed 
RF/FSO systems, atmospheric turbulence may lead to a sig¬ 
nificant degradation of the end-to-end performance. Motivated 
by this limitation, in this paper, we consider a mixed RF 
and hybrid RF/FSO system where an additional RF backhaul 
link is employed as a back-up for the FSO backhaul link. 
Thereby, we assume that the back-up RF link for the relay- 
destination hop utilizes the same bandwidth resource as the 
RF link for the source-relay link. We consider a multi-user 
multiple-input multiple-output (MIMO) setup to fully exploit 
the available RF bandwidth. That is, relay and destination are 
equipped with multiple antennas for data transmission over 
the RF links. For simplicity and feasibility reasons, the relay 
is assumed to operate in the half-duplex mode with respect 
to the RF links, i.e., it can either receive from the users or 
transmit to the destination over the RF links. Therefore, we 
develop an optimal resource allocation policy which allocates 
the RF bandwidth resource to the two RF links based on the 
channel state information (CSI) of the RF and FSO links. The 
proposed protocol adaptively switches between transmission 
and reception for the RF links which requires the relay 
to be equipped with a buffer Q, 0 to temporarily store 
the data received from the users and forward it later over 
the FSO and/or RF backhaul links to the destination. Our 
numerical results confirm the effectiveness of the proposed 
system architecture and resource allocation policy. 

Notations: We use the following notations throughout this 
paper: E{-} denotes expectation, | • | represents the magnitude 
of a complex number and the determinant of a matrix, Z 
denotes the phase of a complex number, and erf(-) is the 
Gauss-error function. Bold capital and small letters are used 
to denote matrices and vectors, respectively. denotes the 
Hermitian transpose of A, I„ is an n x n identity matrix, 
diagjai, 02 ,..., a„} is an n x n diagonal matrix whose 
diagonal elements are 01 , 02 ,... ,a„, and [AJ^n denotes the 
element in the m-th row and n-th column of matrix A. 
!{•} G {0,1} is an indicator function which is equal to one 
if the argument is true and equal to zero if it is not true. 
Additionally, Rice(H,'I'), Unif(a, 6), and GGamma(a,/?) 
denote a Ricean random variable (RV) with parameters H and 
fl', a RV uniformly distributed in the interval [o, b], and a 
Gamma-Gamma RV with parameters a and /3, respectively. 
For notational convenience, we use the definitions C{x) = 
log 2 (l + x), [ccJdB = lOlogiQX, and [a:]'*’ = max{0,a:}. 

II. Preliminaries and Assumptions 

In this section, we describe the considered system model, 
introduce the model for the mixed RF and hybrid RF/FSO 
communication links, and specify our assumptions regarding 
the CSI knowledge. 


signals at the relay and the destination can be modelled as 
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-Radio Frequency (RF) Link - Free Space Optic (FSO) Link 

Fig. 1. Proposed mixed system which consists of cascaded RF and hybrid 
RF/FSO links in a dual-hop configuration. The setup includes K single¬ 
antenna users tlk, k = 1,.. ., K, relay node 7 ?., and destination X>. 

A. System Model 

The considered communication setup is schematically 
shown in Fig. [T] In particular, K users Uk, k = 1,... ,K, 
wish to communicate with destination X> via relay node TZ. 
There is no direct communication link between the users 
and the destination, i.e., the users can send their data to the 
destination only through the relay node. There are two types 
of communication links in our system model: i) the U-TZ 
and TZ-T> RF links and ii) the 7Z-T> FSO link. All RF links 
are assumed to use the same frequency. We assume that each 
user has a single antenna while the relay and the destination 
are equipped with J > K and L antennas, respectively. 
Furthermore, we consider the practical half-duplex constraint 
for the relay node, i.e., the relay cannot receive from the users 
through the U-TZ RF fronthauFaccess link and transmit to the 
destination via the RF backhaul link simultaneously. 

Since the FSO link does not interfere the RF links, the 
relay can always transmit data to the destination via the FSO 
backhaul link. 

Let Tutz, and denote the coherence time of 

the U-TZ RF link, the 7Z-T> RF link, and the 7Z-T> FSO link, 
respectively. Since the users are mobile nodes whereas both 
the relay and the destination are fixed nodes, the U-TZ RF link 
fluctuates considerably faster than the 'TZ-'D RF link and the 
FSO link. Hence, we assume that Tutz ^ holds. 

Thereby, we consider the coherence time of the U-TZ RF link 
as the time reference for resource allocation. In particular, we 
assume the entire time of operation is divided into B blocks of 
length Tutz where each block consists of N symbol intervals 
for the RF signal and each node transmits a codeword which 
spans one block or a fraction of one block. Furthermore, we 
assume that the FSO link has a larger bandwidth than the RF 
links. In order to model this, we assume that the duration of 
the symbol intervals of the FSO signal is M times smaller 
than that of the RF signals, i.e., the FSO signal has a symbol 
rate that is M times larger than the symbol rate of the RF 
signal. Additionally, we assume that users always have enough 
information to send in all blocks and that the number of blocks 
satisfies B ^ oo. 


y-[b]=n,[b]Tc][b]+z][b], j = l ,2 ( 1 ) 

where x"[ 6 ] e and [ 6 ] € denote the transmit 

symbol vectors of the users and the relay, respectively, y" [ 6 ] € 
y 2 [&] G are the received symbol vectors at 

the relay and the destination, respectively, and z"[&] G 
and Z 2 [b] G denote the noise vectors at the relay and 

the destination, respectively, in the n-th symbol interval of 
the &-th block. We assume that z" [&] and Z 2 [&] are zero-mean 
complex additive white Gaussian noise (AWGN) vectors with 
covariance matrices a^Ij and respectively, and are 

mutually independent and independent from the transmitted 
codewords. The variance of the noise at the RF receivers is 
given by [(7|]dB = + Nj^p, j G {TZ,T>}, where 

ppRF, TVj. and denote the RF bandwidth, the noise 
power spectral density (in dB/Hz), and the noise figure of 
the receiver, respectively. Furthermore, Hi[ 6 ] G and 

H 2 [ 6 ] G denote the channel coefficient matrices of the 

U-TZ and TZ-'D RF links, respectively. Moreover, we assume 
that all elements of Hi [ 6 ] and H 2 [ 6 ] are mutually independent, 
ergodic, and stationary random processes with continuous 
probability density functions. Throughout the paper, we as¬ 
sume the RF links are affected by Ricean fading 111, IH. We 
further take into account the effect of path loss. Suppose h[b] 
is one of the elements of Hj[5], j = 1,2. Thereby, h[b] can 
be written as h[b] = \/h^hf[b], where ha and hf[b] are the 
average power gain and the fading coefficient of the RF link 
which are given by a, 13 


ha — 




47rd« 


“ref 


( 2 ) 


|/i/[&]| ^ Rice(H,T'), Zhf[b] ~ Unif(— tt,tt) 


where is the wavelength of the RF signal, and G^^ 
are the RF transmit and receive antenna gains, respectively, 
is a reference distance for the antenna far-held, is 
the distance between RF transmitter and receiver, and z/ is 
the RF path-loss exponent. Parameters H and in the Rice 
distribution are the ratio between the power in the direct path 
and the power in the scattered paths and the total power in both 
paths, respectively. For the U-TZ RF link, we assume that the 
direct link is not available (H = 0). Hence, the Ricean fading 
reduces to Rayleigh fading. Finally, the transmitted RF signals 
have to meet the following power constraints 


E{xn 6 ](x?[ 6 ])^} < dmg{P^,P^, ...,P^) (3a) 

E{(x-[&])^x-[ 6 ]} < (3b) 


where P^ and P^^ are the RF transmit powers of user k and 
the relay, respectively. 


B. Communication Links 

In the following, we describe the channel models for the 
RF and FSO links which are used throughout the paper. 

RF Links: The U-TZ RF link constitutes a {K, J)-distributed 
MIMO system with K single-antenna transmitters and one J- 
antenna receiver. In contrast, the TZ-D RF backhaul link is 
a standard (J, P)-point-to-point MIMO system. The received 


FSO Link: The relay node is equipped with an aperture 
transmitter pointing in the direction of a photo detector at 
the destination. We assume an intensity modulation direct 
detection (IM/DD) FSO system with on-off keying (OOK) 
modulation for the FSO link, i.e., the photo detector directly 
detects the intensity of the received photo-current by integrat¬ 
ing over each symbol interval. In particular, after removing 
the ambient background light intensity, the detected signal 


















A. The Proposed Protocol 


intensity at the destination is modelled as 0 


y^[b]=pg[b]x^[b]+z"^[b], (4) 


where p is the responsivity of the photo detector and a;™ [ 6 ] € 
{O, G R, and G R are the power 

of the OOK modulated symbol at the relay, the detected 
signal power at the destination, and the shot noise at the 
destination caused by ambient light for the m-th symbol of 
the 6 -th block, respectively. Moreover, is the maximum 

allowable transmit power of the relay over the FSO link which 
is mainly determined by eye safety regulations H]. Noise 
z™[ 6 ] is modelled as zero-mean real AWGN with variance 
and is independent of the transmitted signal. Moreover, 
g[b] is modelled as g[b] = 5 aff/[ 6 ], where ga and gf[b] are the 
average gain and the fading gain of the FSO link, respectively, 
and are given by 



-y/^r 

GGamma(a, /3) 


X 10 


— Kd^ 


^/lO 


(5) 


where r is the aperture radius, (f> is the divergence angle 
of the beam, is the distance between FSO transmitter 

and receiver, and k is the weather-dependent attenuation 
coefficient. Assuming spherical wave propagation, parameters 
a and /3 in the Gamma-Gamma distribution are related to 
physical parameters as follows a 


a = 


exp 


/3 


exp 


0.49t?2 

[1-f 0.18^2-f0.56t?i2/5]^/6 


-1 


- 1 


( 6 a) 


0.51i? 2 [l + 0.69t?12/5j 5/6 I 

(1 + 0 . 9^2 + 0 . 62 ^ 2 ^ 12 / 5 ) 5/6 J 


, ( 6 b) 


where ^ and 

C = 27r/A^®°. Here, is the wavelength and is the 

weather-dependent index of refraction structure parameter. 


C. CSl Knowledge 

Throughout this paper, we assume that the relay has full 
knowledge of the CSI of all links and is responsible for de¬ 
termining the optimal transmission strategy and for conveying 
the strategy to the other nodes, cf. Theorem[T] Moreover, since 
for large K, an excessive amount of CSI feedback is required 
for adaptive rate transmission at the user nodes, we assume 
that the users transmit with a priori fixed transmission rates 
and hence, no CSI knowledge is required at the users. We 
note, however, that the proposed problem formulation and the 
resulting protocol are provided in a general form and the case 
where the users transmit with adaptive rates can be easily 
accommodated, cf. Remark [T] Furthermore, the destination 
knows the CSI of the 'R.-'D links as required for reliable co¬ 
herent decoding. Moreover, we assume that the channel states 
change slowly enough such that the signaling overhead caused 
by channel estimation and feedback is negligible compared to 
the amount of information transmitted in one block. 


III. Mixed RF and Hybrid RF/FSO Protocol 
In this section, we first introduce the considered protocol 
and explain the adopted coding scheme. Subsequently, we 
derive the optimal resource allocation policy for the proposed 
protocol. 


Recall that due to the half-duplex constraint, the relay 
cannot use the U-TR and 'R.-'D RF links simultaneously. In 
fact, if the quality of the FSO backhaul link is sufficiently 
good, all data received at the relay from the users can be 
forwarded using the FSO link to the destination, and there is 
no need to activate the 'R.-D RF backhaul link at all. However, 
due to the time-varying fluctuations of the FSO link, the relay 
may not always be able to forward the users’ information to 
the destination by employing only the FSO link. Hence, the 
7?.-X> RF link is needed as a back-up. On the other hand, using 
the 'R.-D RF backhaul link comes at the expense of reducing 
the time of transmission of the users to the relay over the lA-'R 
RF link due to the half-duplex constraint. 

In light of the above discussion, the main idea of the 
proposed protocol is to adaptively share the RF bandwidth 
resource between the lA-'R RF fronthaul and the R-D RF 
backhaul link based on the fading state in the 6 -th block, 
i.e.. Hi[ 6 ], H 2 [ 6 ], and g[b]. To this end, we introduce time 
sharing variable q\b] G [0,1] where q\b] denotes the fraction 
of block 6 in which the relay receives, i.e., the U-R RF link is 
active. In the remaining 1 — ( 7 [ 6 ] fraction of block 6 , the relay 
transmits over the R-D RF link. We assume that the relay 
node is equipped with an infinite-size buffer for data storage. 
Let (5[6] denote the number of information bits available in 
the buffer after receiving the information from the users in 
the 6 -th block. Using these notations, the coding schemes, 
the transmission rates, and the dynamics of the queues at the 
buffers are presented below. 

lA-R RF Link: The users employ Gaussian codebooks, i.e., 
the fc-th element of vector Xi[ 6 ] is generated independently 
according to a zero-mean rotationally invariant complex Gaus¬ 
sian distribution with variance P^. At the beginning of each 
fading block, each user k encodes q[b]NR^ bits of information 
into a codeword taken from a Gaussian codebook with a fixed 
rate bits/symbol and having a length of q[b]N symbols. 
The users transmit their codewords and the relay receives 
y" [ 6 ] according to ([U. The relay can employ several multi¬ 
user detection schemes proposed in the literature El, e.g., 
linear zero-forcing (ZF) and minimum mean square error 
(MMSE) detection, nonlinear detection schemes incorporating 
successive interference cancellation (SIC), or optimal maxi¬ 
mum likelihood (ML) detection. Depending on the type of 
detector used, the codeword of each user experiences a certain 
signal-to-noise ratio (SNR) denoted by 7 ^ [ 6 ]. For instance, for 
ZF detection, 7 ^ [ 6 ] is given by 


7£[&] = 


rtk 


"'R. 


((Hi[6])«Hi[6]) 


-1 


(7) 


- kk 


The codeword of user k can be decoded reliably only if 
> 2 ^ — 1 ; otherwise the relay cannot decode this 
codeword and has to drop it and ask user k to retransmit 
this information in the following blocks. The effective rate 
of the lA-R RF link is defined as the normalized sum 
of the information bits in bits/symbol that can be decoded 
successfully at the relay and is given by 


c^^[b] = J2H'rh{b]>2^^-i}Ru- ( 8 ) 














IZ-T) RF Link: This backhaul link is a standard point-to- 
point MIMO channel with an average power constraint across 
all antennas. We employ Gaussian codebooks, i.e., vector X 2 [t»] 
is a multivariate zero-mean rotationally invariant complex 
Gaussian vector, and waterfilling power allocation across the 
transmit antennas Go). We note that the relay cannot transmit 
more information than what it has stored in its buffer and has 
not yet been transmitted over the FSO link. Thereby, the relay 
extracts min{Q[5],[l — g[ 6 ]]A^C^^[ 6 ]} bits of information 
from its buffer and transmits them over its J antennas where 


C. 


RF 


min{ J,L} 

[‘I" E 

i=i 


log 2 M 


4 

.2 

T> 


T + 


(9) 


Here, Xj is the j-th singular value of H 2 ifTOl : /i is the water 
level which is chosen to satisfy the power constraint in (l^^) 
as the solution of the following equation 


min{ J,L} 




^2? 

x]j 


_ pRF 


( 10 ) 


FSO Link: We employ OOK modulation and soft detection 
for the FSO link. From an information theoretical point of 
view, the FSO link can be modelled as a binary input- 
continuous output AWGN channel where the capacity is 
achieved with a uniform distribution of the binary inputs im, 
im. Thereby, at the beginning of each block, the relay extracts 
min { [Q[b- 1] - (1 - q[b- l])iVC' 2 ^i^[ 6 - 1]] + , MNC^^°[b]} 
bits of information from its buffer, encodes them into a 
codeword, and sends the codeword to the destination where 
(^FSOj^j is the capacity of the Gaussian channel with OOK 
input and is given by im, ina 


C^^^[b] = 1 - ^ exp {-f 2 } log 2 I 


1 -I- exp I 

_^v 

2ct2 / 


2ct2 / 

jdf, (11) 


where p[b] = pg[b]P^^. 

Dynamics of the Queue: After the relay has received the 
information from the users in the &-th block, the amount of 
information bits in the buffer is updated as 


Finding the maximum throughput and the corresponding op¬ 
timal resource allocation policy is challenging for the con¬ 
sidered channel model due to the recursive dynamics of the 
queue. Nevertheless, we can obtain an upper bound on the 
achievable throughput of the proposed protocol by neglecting 
the effect of the queues on the transmission rate in (fTSI) . 

Upper Bound The achievable average rate of the 

proposed protocol for the considered mixed RF and hybrid 
RF/FSO system with block fading is upper bounded by t'^pp 
obtained from the following optimization problem 

maximize t'^pp subject to (14) 

<j'[&]e[o,i], Vf) 

1 ° 

Cl : r"PP < lim — V g[ 6 ]lVCP[ 6 ] 

1 ® 

C2 : r'^PP < ^ [(1 - q[b])NCf^[b] + NMC^^^[b]]. 


In the following theorem, we provide the optimal solution 
to the above optimization problem. For notational simplicity, 
let C'P^^(Hi), C'^^(H 2 ), and denote the capacity 

functions given in ([ 8 ]l, (|9]l, and (fTTl i in terms of the fading 
states Hi, H 2 , and g, respectively. Moreover, let /hi(Hi), 
/h 2 (H 2 ), and fg{g) denote the probability density functions 
of Hi, H 2 , and g, respectively. 

Theorem 1: The optimal resource allocation policy depends 
only on the values of channel matrices Hi and H 2 and 
is denoted by ( 3 '*(Hi,H 2 ) G [0,1]. Moreover, the optimal 
resource allocation policy reduces to the adaptive activation 
of either the 14-71. link or the 7?.-X> RF link according to 


9*(Hi,H2) 


0, ii 

1 , otherwise 


where A G (0,1] is a constant which does not depend on 
the instantaneous realization of the fading but depends on the 
fading distributions. The optimal value of A can be obtained 
offline before transmission starts using an iterative algorithm 
with the following update equation 


A[* -I- 1] 






1 

0 


, (16) 


Q[b] 


Q[b-l]-il-q[b-l])NCf^[b-l] 




q[b]NCf^[b]. (12) 


B. Optimal Resource Allocation 


In this section our goal is to obtain the optimal resource 
allocation for the RF links, i.e., the optimal q[b], V 6 , such 
that the average number of information bits received at the 
destination, denoted by r, is maximized. Note that since there 
is no direct link, the amount of information received from the 
users at the destination is identical to that received from the 
relay at the destination. Hence, the throughput maximization 
problem can be formulated as 


where i is the iteration index and (5[z] is an appropriately 
chosen step size parameter. The average capacity rates 
and (7^®° are given by 

= E{q*(Hi,H2)Cr(Hi)} 

= JJ g*(Hi,H2)Cr(Hi) 

x/Hi(Hi)/H 2 (H 2 )dHidH 2 (17a) 

= E{[1 - r7*(Hi,H2)]C2^^(H2)} 

= |y'[l-g*(Hi,H2)]C'2^^[H2] 

x/Hi(Hi)/H2(H2)dHidH2 (17b) 

^ = I C^^°{g)fg{g)dg, (17c) 


’ ’ b—1 

-fmin {[Q[b-l] - {l-q[b-l])NCf^[b-l]]+ [&]}. 


where for a given A = A[*] in the *-th iteration, g*(Hi, H 2 ) has 
to be obtained from (fTSl l for a given set of fading values. Em¬ 
ploying the optimal A* obtained from the iterative algorithm 
and the optimal q* (Hi, H 2 ), from (fTsT l. and substituting them 

















into (O, the upper bound is obtained as 

r'^PP = iVmin{Cf^,C^^( 18 ) 

Proof: Please refer to Appendix lAl ■ 

In the following, we formally state the achievability of the 
upper bound t'^pp as B ^ oo. 

Lemma 1: The upper bound t“pp is achievable if the op¬ 
timal resource allocation policy in Theorem [T] is employed. 
More precisely, as B ^ oo, we obtain r —>■ t'^pp. 

Proof: Please refer to Appendix iBl ■ 

We note that Lemma [1] reveals that t'^pp is asymptotically 
achievable by using the resource allocation policy in Theo¬ 
rem [T] In fact, it is possible that for some fading blocks the 
relay does not have enough information in its queue to send 
to the destination over the FSO and 7?.-X> RF links, but as 
B —> oo, the effect of these events on the throughput becomes 
negligible, i.e., r —^ r“PP. 

Corollary 1: The optimal value of Lagrange multiplier A 
can be obtained analytically from one of the following mutu¬ 
ally exclusive cases: 

Case 1: If E{C'f^^(Hi)} < holds, we obtain 

A* = 1. Hence, we have g*(Hi,H 2 ) = 1, VHi,H 2 . 

Case 2: If E{C'i^P(Hi)} > ME{C^®°(p)} holds, the 
optimal value of A is obtained analytically from the following 
equation 

Jj <7*(Hi,H2)Cr(Hi)/Hi(Hi)/H,(H2)dHidH2 
= yy[l-g*(Hi,H2)]Ci^^(H2)/Hi(Hi)/H,(H2)dHidH2 

+M I C^®°(5)/g(5)d5, (19) 

where ( 3 '*(Hi,H 2 ) is substituted from (flSl l. 

Proof: Please refer to Appendix ICl ■ 

Corollary [T] provides an important intuition based on the 
statistics of the RF and FSO links. In particular, if the FSO 
link is statistically strong (cf. Case 1), the relay can forward all 
the information received from the users to the destination over 
the FSO backhaul link and the 'TL-'D RF backhaul link remains 
inactive for all fading blocks, i.e., we obtain A* = 1. However, 
if the statistical quality of the FSO link is not as strong e.g. due 
to severe atmospheric conditions (cf. Case 2), the 'TL-'D RF 
backhaul link becomes active in fading blocks which satisfy 
ACf^^(Hi) < [1 — A]C'|^^(H 2 ). Thereby, the optimal value 
of A is obtained such that the sum of the average information 
rates sent from the relay to the destination over both the RF 
and the FSO links is identical to the average information rate 
received at the relay from the users. 

Remark 1: We note that the proposed problem formulation 
and the resulting protocol in Theorem [T] are provided in a 
general form such that coding schemes different from those 
considered in this paper can be easily accommodated. In 
particular, using different coding schemes changes only the 
values of the instantaneous capacity expressions in ®, (|9]l, 
and (HB such that the analysis of the proposed protocol 
and resource allocation policy remains valid. For instance, as 
explained in Section II.C., we assume that the users transmit 
with a priori fixed transmission rate in order to avoid the CSI 
feedback overhead from the relay to all the users. However, if 
this CSI overhead can be accommodated, then the users can 


TABLE I 

Values of the Numerical Parameters (2), O- 


RF Link 

Symbol 

Delinilion 

Value 


Distance between user k and the relay 

400 m 

Pu 

Transmit power of user k 

20U mW (23 dBm) 


Relay transmit power over RF link 

2 W (33 dBm) 

d 

Disiiince between the relay and the destination 

|1, 2| km 

d^f, 

Reference distance of the RF link 

60 m 


Antenna gains for U-'R. RF link 

Antenna gains for TZ-T) RF link 

(0,8) dBi 
(10,15) dBi 

No 

Noise power spectral density at the RF receivers 

-114 dBm/MHz 

Np 

Receiver noise ligure of the RF receivers 

5 dB 


Wavelength of RF signal 

85.7 mm (3.5 GHz) 

W’^ 

Bandwidth of RF signal 

20 MHz 

(n,®) 

Rice distribution partimeters for U-TZ RF link 

(0,1) 


Rice distribution parameters for TZ-T) RF link 

(4,1) 

u 

RF path-loss exponent 

3.5 

FSO Link 

Symbol 

Definition 

Value 

d 

Distiince between the relay and the destination 

[1,2J km 


Relay transmit power over FSO link 

40 mW (16 dBm) 


Noise variance at the FSO receivers 

10-^* 


Wavelength of FSO signal 

1550 nm (193 THz) 


Bandwidth of RF signal 

1 GHz 

P 

Responsivity of FSO photo-detector 


K 

Weather-dependent attenuation factor for 

10-4 


[clear air. haze, light foe. moderate foe. heavv foel 

[0.43,4.2,20,42.2,125] 


Weather-dependent index of refraction structure for 
[clear air. haze, light fog, moderate fog, heavy fog] 

10-1® X [50,17,3,2,1] 

4> 

Laser divergence angle 

2 mrad 

r 

Aperture radius 

10 cm 


transmit with adaptive rates in each fading block. Thereby, 
the maximum sum capacity of the multiple-access channel is 
given by E] 


=log2 


Ij + 


1 


xHi[6]diag{Pi,P^,...,P2^}(Hi[6]) 


\H 


■ (20) 


IV. Numerical Results 


The values of the parameters for the RF and FSO links 
used for the numerical results shown in this section are 
provided in Table I. Moreover, we consider modified versions 
of the following two protocols as benchmark schemes i) mixed 
RF/FSO relaying which does not employ the RF backhaul 
link ||6l and ii) conventional RF relaying ifTJI which does not 
employ an FSO backhaul link. Unlike the original versions 
of these protocols, the relays of the benchmark schemes are 
also equipped with buffers. Hence, the nodes are able to 
transmit data for many consecutive fading blocks to average 
out the fading. Additionally, for conventional RF relaying, we 
optimize the ratio of the number of fading blocks allocated 
to the lA-lZ and TZ-'D RF links for throughput maximization, 
similar to Q. Furthermore, for all considered protocols, we 
assume that the relay employs the ZF detector to recover the 
information sent by the users. 

In Fig. m we show the normalized achievable sum through¬ 
put of all users f = (in bits/second) versus the weather- 
dependent attenuation coefficient of the FSO link, k, for 
AT = 5, J = L = 10, d = 1 km and 2 km, and = 8 
bits/symbol for Vfc. From Fig. |2l we observe that for favorable 
atmospheric conditions, i.e., low values of k, the relay is able 
to forward all the information received from the users to the 
destination over the FSO link. Hence, the 'TZ-'D RF backhaul 
link remains inactive in the proposed protocol, cf. Case 1 in 
Corollary [T] As a result, the achievable throughput of the 
proposed protocol and the mixed RF/FSO relaying protocol 
coincide. However, as k increases, i.e., the quality of the FSO 
link deteriorates due to severe atmospheric conditions, there 
is a critical value of k above which the FSO link is no longer 
able to forward all the user data, which is received at the 
relay, to the destination. Hence, the achievable throughput of 
the mixed RF/FSO relaying protocol decreases ultimately to 
zero for large values of n. In contrast, although the achievable 



































Weather-dependent attenuation coefficient of the FSO link, k, (xlO 

Fig. 2. Normalized throughput, f, in Mbits/second vs. weather-dependent 
attenuation coefficient of the FSO link, k, for K = b, J = L = 10, d = 1 km 
and 2 km, and = 8 bits/symbol Vk. The remaining parameters are given 
in Table I. Moreover, from low to high values of k, the vertical dashed-dotted 
lines represent the following weather conditions O: clear air, haze, light fog, 
and moderate fog, respectively. 


multipliers corresponding to constraints Cl and C2 by Ai and 
A 2 , respectively, the Lagrangian function corresponding to the 
optimization problem in (fT4l t is obtained as 


£((?[&], r"PP,Ai,A2) =t"pp 


+Ai 




b=l 


( 21 ) 


-I-A 2 


^ ^ [[1 - q[b]]NCf^[b] + MNC^^°[b]] - r"PP 


The dual function is then given by 


r’(Ai,A 2 )= maximize £((;[6], t'^pp, Ai, A 2 ) (22) 

i3[f)]G[0,l],r“PP>0 

and the dual problem is given by 


minimize 22(Ai, A 2 ). (23) 

Ai>0,A2>0 

To solve (fT4l i using the dual problem in (l23b . we first obtain 
the primal variables q[b] and r“PP for a given dual variables 
Ai and A 2 from (l22l i. Then, the optimal Ai and A 2 are obtained 
by solving the dual problem in (|2^ . 


throughput of the proposed protocol also decreases due to the 
sharing of the RF bandwidth between the lA-Tt and 7t-T> RF 
links, it converges to a non-zero value for large values of k, cf. 
Case 2 in Corollary [T] Note that the achievable throughput of 
the conventional RF relaying protocol does not depend on the 
quality of the FSO link. Moreover, for high values of k when 
the FSO link becomes unavailable, the proposed protocol still 
outperforms the conventional RF relaying protocol due to the 
adaptive resource allocation based on the instantaneous CSI, 
cf. Theorem [U Finally, we observe that the critical value of k 
above which adaptive RF link selection becomes necessary 
depends on the system parameters. For instance, for the 
parameters considered here and moderate foggy atmospheric 
conditions, i.e., k = 20 x 10“^, the 'TL-'D RF link is active if 
d = 2 km but inactive if d = 1 km. 

V. Conclusions 

A mixed RF and hybrid RF/FSO system was considered 
where multiple users transmit their data over an RF link to 
a relay node and the relay forwards the information to the 
destination over a hybrid RF/FSO backhaul link. The optimal 
resource allocation policy which shares the RF bandwidth 
resource between the RF links based on the CSI was de¬ 
rived. Our numerical results revealed the effectiveness of the 
proposed communication architecture and resource allocation 
policy even if the FSO link is affected by severe atmospheric 
conditions. 


Appendix A 
Proof of Theorem[T| 

In this appendix, our aim is to find the optimal resource 
allocation policy as a solution of the optimization problem 
given in Gil. Since the cost function and the constraints in 
d are affline in the optimization variables q[b], Mb, and 
and the feasible set is non-empty. Slater’s condition is satisfied. 
Hence, the duality gap is zero Gl- Therefore, the solution of 
the primal problem in Gl can be found from the solution 
of the dual problem of d d- Denoting the Lagrange 


A. Optimal Primal Variables 

The optimal resource allocation variables, q\b], Mb, and 
the optimal value of the upper bound, are either the 

boundary points of the feasible sets, i.e., q]})] € [0,1] and 
r'^PP > 0, or the stationary points which can be obtained 
by setting the derivatives of the Lagrangian function in (I 2 TI 1 
with respect to r“PP and to zero. The derivatives of the 
Lagrangian function in (l2TT i are given by 


_ 1 _ Ai - A 2 

LL = I [A.jvcrii] - A,«cn(,]] . 


(24a) 

(24b) 


If the derivative is non-zero, the optimal value of t“pp 

is at the boundary of its feasible set, i.e., t'^pp —>■ 00 or 
t“pp —0, which cannot be the optimal solution. Therefore, 
the derivative in d^l l has to be zero which leads to 

Ai + A 2 = 1. Since the probability density functions of the 
channel coefficients Hi and H 2 are continuous, we obtain that 
Pr = o| = 0 holds. Hence, the optimal value of q\b] is 

always at the boundaries, i.e.. 




1 , if XiCf^[b]>X 2 Cf^[b] 

0, otherwise 


(25) 


Defining A = Ai = 1 — A 2 , (l25l l leads to the optimal link 
selection policy in G3]) . 


B. Optimal Dual Variable 

The optimal value of A is obtained by solving the dual 
problem in (l23t . In particular, substituting, the optimal values 
of q*[b] in ( l25b for a given A into (I 2 TI 1 . we obtain the dual 
function in (l22l i. The dual function D{X) is a linear function 
of A and its derivative is obtained as 


6=1 


-[(1 - -f MNC^^‘^[b]] 


■(26) 








































We note that for a given A, the value of is fixed. 

Thereby, we can investigate the following mutually exclusive 
cases for the optimal A*; 

Case 1: If we assume that A* > 1 holds, we obtain 
= 1, V6 from (l25l l. Thereby, the derivative becomes 
positive and in order to minimize 'D{\), we have to decrease 
A. In fact, for each A* > 1 and e > 0, we can find a new 
A = A* —e where q]})] = 1, V6 still holds and the dual function 
is decreased. Hence, A* > 1 cannot be optimal. 

Case 2: If we assume that A* = 0 holds, we obtain q[b] = 
0, \/b from (l25T l. Intuitively, this cannot be the optimal dual 
variable since it leads to = 0. In particular, the derivative 
becomes negative and in order to minimize T’(A), we 
have to increase A. In fact, we can assume a small A = e and 
obtain a positive t“pp and decrease the dual function. Hence, 
A* = 0 cannot hold for the optimal solution. 

Case 3: Considering Case 1 and Case 2, 0 < A* < 1 has 
to hold. In order to find the optimal dual variable, we can 
employ the widely-adopted gradient method Ea. The main 
idea is to minimize 22(A) by updating A along the gradient 
search direction. The updates are performed as 


A[i + 1] = 


A[i] — 5[i 


522(A) 


dX 


A=A[2] 


(27) 


This leads to the iterative algorithm in Theorem [T] with the 
updates given in (fThl) . The gradient method is guaranteed to 
converge to the optimal dual variable A* provided that the step 
sizes ()[i] are chosen sufficiently small ifTsl . This completes the 
proof. 


Appendix B 
Proof of Lemma[T] 

The proof follows from a well-known property in queuing 
theory M- Let us consider a buffer with stochastic arrival 
process a[i], departure demand process (queue process) d\i], 
and queue lengths Q[i\ in the i-th time instance. Thereby, 
the departure rate (queue service) is obtained as b[i] = 
min{Q[i — l],(2[i]} and Q[i] evolves as Q[i] = max{Q[i — 
1] — d[i], 0} -I- a[i]. However, the average departure rate E{&} 
(averaged over the time slots) can be expressed independently 
from the dynamics of the queue as 

{ E{d}, if E{a} > E{d} 

E{a} = E{d}, if E{a} == E{d} (28) 
E{a}, if E{a} < E{d} 

In particular, if E{a} > E{(i} holds, i.e., the average rate 
flowing into the buffer is larger than the average capacity 
of the respective departure channel, then, as the number of 
time instances grow to infinity, the buffer always has enough 
information to supply because the amount of information in 
the queue increases over time and we obtain E{&} = E{(2}. 
On the other hand, if E{a} < E{(2} holds, i.e., the average 
information flowing into the buffer is less than the average 
capacity of the respective departure channel, then by the 
law of conservation of flow, we obtain E{6} = E{a}. If 
E{a} = E{(i} holds, we obtain E{6} = E{a} = Ejd} due to 
the continuity property. Hence, from (l28l) . we can conclude 
that E{6} = min{E{a}, E{d}} holds which leads to the 


expression given in (flSl l for the upper bound, i.e., r = 

This completes the proof. 

Appendix C 

Proof of Corollary [T] 

If E{C'f’^(Hi)} < ME{C'P®0(g)} holds, from we 
obtain that ^ < 0 holds regardless of the value of A. 
Hence, in order to minimize the dual function 22(A), we 
have to increase A. Moreover, since A € (0,1] holds, we 
obtain A* = 1. On the other hand, if E{C]^^(Hi)} > 
ME{C'f®°(p)} holds, we can conclude that 0 < A* < 1 and 
its value is unique. In other words, according to (fTSl l. = 
N min where is a monotonically 

increasing function of A, C 2 is ^ monotonically decreasing 
function of A, and does not depend on A. Therefore, the 
optimal A has to be chosen such that 
holds which leads to (O. This completes the proof. 
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